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We trace the development of azimuthal anisotropy n = 2,3) via parton-parton collision history 
in two transport models. The parton is studied as a function of the number of collisions of each 
parton in Au+Au and d+Au collisions at = 200 GeV. It is found that the majority of 

comes from the anisotropic escape probability of partons, with no fundamental difference at low 
and high transverse momenta. The contribution to Vn from hydrodynamic-type collective flow 
is found to be small. Only when the parton-parton cross-section is set unrealistically large does 
this contribution start to take over. Our findings challenge the current paradigm emerged from 
hydrodynamic comparisons to anisotropy data. 
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Relativistic heavy ion collisions aim to create the 
quark-gluon plasma (QGP) and study quantum chromo¬ 
dynamics (QCD) at the extreme conditions of high tem¬ 
perature and energy density [1]. The system created in 
these collisions is described well by hydrodynamics where 
the high pressure buildup drives the system to expand 
at relativistic speed [5]. Experimental data fit with hy¬ 
drodynamics inspired models suggest that particles are 
locally thermalized and possess a common radial flow 
velocity [3]. Of particular interest are non-central col¬ 
lisions where the overlap volume of the colliding nuclei 
is anisotropic in the transverse plane (perpendicular to 
beam). The pressure gradient would generate anisotropic 
expansion and final-state elliptic flow [3]. Large ellip¬ 
tic anisotropy in momentum (U 2 ) has been measured, as 
large as hydrodynamic calculations predict [U |2] . This 
suggests that the collision system is strongly interacting 
and nearly thermalized (sQGP) j^. 

Molnar’s Parton Gascade (MPG) [IT] can describe the 
measured large V 2 , but with an unusually large parton- 
parton interaction cross-section (cr) [S]. It approaches [7] 
the limiting case of ideal hydrodynamics (cr —> 00 ) and 
may be an effective description of the sQGP. A Multi- 
Phase Transport (AMPT) [51 [S] can describe the large 
anisotropy with tr motivated by perturbative QGD but 
with the string melting mechanism m- String melting 
liberates strings into a larger number of quarks and an- 
tiquarks, effectively increasing the system opacity m- 

It is generally perceived that large V 2 can only be gen¬ 
erated in large-system heavy ion collisions. Recent parti¬ 
cle correlation data, however, hint at similar V 2 in small 
systems of high multiplicity p+p and p-|-Pb collisions at 
the LHC [Hj and d-|-Au collisions at RHIC [T3]. Hydro¬ 
dynamics has been applied to these systems and seems 
to describe the experimental data well [14] . AMPT also 
appears to describe the measured correlations [15j . This 


suggests that these small-system collisions might create 
an sQGP as well, in contrast to the general expectations. 

The purpose of this Letter is to study the development 
of azimuthal anisotropy to shed light on its connection 
to the properties of the sQGP and thermalization. We 
employ the string melting version of AMPT [9] because it 
reasonably reproduces particle yields, transverse momen¬ 
tum (p_l) spectra, and V 2 data for the bulk (see Figs. 1-3 
of Ref. ffH). The model consists of a fluctuating initial 
condition, parton elastic scatterings, quark coalescence 
for hadronization, and hadronic interactions. In partic¬ 
ular, parton scatterings are treated with Zhang’s Parton 
Gascade (ZPG) [15]. We use Debye screened differential 
cross-section da/dt oc a1/{t — [5] in AMPT, with 

strong coupling constant as = 0.33 and Debye screen¬ 
ing mass pd = 2.265/fm (the total cross section is then 
CT = 3 mb). We also employ MPG [17], in order to check 
the generality and model dependence of our study. MPG 
employs the parton subdivision technique to eliminate 
acausal numerical artifacts due to action at a distance in 
the cascade algorithm. Here MPG is used with smooth, 
longitudinal boost invariant, binary collision profile for 
the initial conditions, as in Ref. [7]. We use both the 
Debye and isotropic da/dt in MPG; the results are quali¬ 
tatively similar. We analyze the entire history of parton- 
parton interactions in these models. For simplicity only 
partons are analyzed. 

We compute the harmonic plane (short-axis direc¬ 
tion of the corresponding harmonic component) of each 
event from its initial configuration of partons |18j by 


= — \dXd,n2{{r\ sinnc/r) 1 {r\ cosnc/r)) + tt] . (1) 


Here rj^ and ipr are the polar coordinate of each initial 
parton (after its formation time) in the transverse plane, 
and (...) denotes the per-event average. We analyze the 


2 



FIG. 1: (a) Normalized probability distributions of partons 
freezing out after Ncoii collisions in AMPT (Debye a = 3 mb). 
The thin and thick solid curves are for partons of all p±_ in 
d+Au and Au+Au collisions, respectively. The dashed curve 
is for partons with final px>3 GeV/c in Au+Au collisions, 
(b) Normalized probability distributions of parton transverse 
radius rx in Au+Au collisions from AMPT. The thick curves 
are for A^coii=0 and the thin curves for A^coii=5. The solid 
curves are for freezeout partons and the dashed curves for 
active partons. 


momentum anisotropy in the initial state, final state, 
and any intermediate state in-between. The momentum 
anisotropy is characterized by Fourier coefficients |19j 

= (cosn(0-, (2) 

where (j) is the azimuthal angle of the parton momentum. 

Results. We simulate Au+Au (impact parameter b = 
7.3 fm) and d+Au (5 = 0 fm) collisions by AMPT with 
CT = 3 mb. We trace the history of parton cascading 
by the number of collisions (A^coii) a parton suffers with 
other partons. Figure [^a) shows the probability distri¬ 
butions of partons freezing out after Ncoii collisions. Par- 
tons in mid-central Au+Au suffer more collisions than in 
d+Au, as expected. See Table |T] for the average number 
of collisions (opacity), (A^cou), for each collision system. 
As seen in Fig. [^a), partons with higher final p± have 
fewer collisions and freeze out earlier. 

As shown in Fig. [^a), some partons do not interact 
at all and thus instantly freeze out at A^coii=0. These 
partons tend to reside in the outer region of the over¬ 
lap volume (“surface emission”), as shown in Fig. [^b) 
for Au+Au where the transverse radius (rx) distribu¬ 
tion of freezeout partons is depicted by the thick solid 
curve. Those continuing to interact tend to be inside as 
shown by the thick dashed curve. This feature is qualita¬ 
tively similar for all A^coii values (e.g. see the thin curves 
in Fig. [^b) for A^coU=5). This is consistent with the 
general expectation-the energy density is smaller in the 
outer shell thus the probability for further interactions 
is smaller. It is interesting to note that the freezeout 
“surface” moves inward, indicating an outside-to-inside 
freezeout scenario. 


FIG. 2: (a) Gomparison of the normalized probability distri¬ 
butions of freezeout partons in Au+Au collisions from AMPT 
(Debye cr = 3 mb) and MPC (both Debye and isotropic 
a = 5.5 mb) of comparable opacities, (b) Normalized proba¬ 
bility distributions of partons freezing out after Acoii collisions 
in MPG with three isotropic a values: 5.5 mb (solid), 20 mb 
(dashed), and 40 mb (dash-dotted). 


Due to different initial parton densities, we use De¬ 
bye a = 5.5 mb for MPC to obtain a similar opacity 
(A);oii)=4-5 as AMPT with Debye cr = 3 mb. The Ncoii 
distributions are similar between AMPT and MPC, as 
illustrated in Fig.j^a). Also shown is the MPC A^coii dis¬ 
tribution with isotropic cr = 5.5 mb which is again simi¬ 
lar. Figure [^b) shows the Woii distributions from MPC 
with three isotropic a values: 5.5 mb, 20 mb, and 40 mb. 
As cr increases, {Ncoii) becomes larger as expected. The 
probability for small Woii is, nevertheless, non-zero even 
at large cr because of finite surface emission. 

We track the development of V 2 in AMPT and MPC 
by studying parton V 2 as a function of Woii in Fig. 
for Au+Au collisions. The solid curves are the V 2 of all 
partons after suffering A^coii collisions (those frozen out 
with smaller A^joii values are of course not included). The 
dashed curves are the V 2 of the partons that freeze out 
after suffering exactly A^coii collisions (i.e. without fur¬ 
ther interactions). At A"coU = 0, the V 2 of all partons 
is zero because parton form with axially uniform mo¬ 
menta. Some partons (“corona”) do not interact at all 
and instantly freeze out with Ncoii = 0. Because there is 
a larger probability for the partons to escape along the 
short axis of the overlap volume, those freezeout partons 
have positive V 2 muin]. In the low density limit (LDL), 
the anisotropy may be analytically derived EH 122]. In 
fact, this escape mechanism is rather general as it hap¬ 
pens throughout the entire evolution of the collision sys¬ 
tem. After Ncoii collisions, the V 2 of all partons is still 
roughly zero. Some of these partons freeze out; they have 
positive V 2 partly due to the preferential escape along the 
short axis. The remaining partons, that can have nega¬ 
tive V 2 as shown by the dotted curves in Fig. [^a), con¬ 
tinue to interact. With one more collision, the azimuthal 
distribution of those partons becomes roughly isotropic 








3 



FIG. 3: Parton V 2 in Au+Au collisions after suffering A^coii 
collisions in (a) AMPT with Debye ct = 3 mb (thick curves) 
and MPC with isotropic a = 5.5 mb (thin curves), and (b) 
MPC with isotropic a = 20 mb. The dashed curves are for 
partons freezing out after exactly Acoii collisions, the dotted 
curves for partons continuing to interact, and the solid curves 
for all partons (i.e. sum of the former two). 


again, with approximately zero V 2 (solid curves). This 
process then repeats itself. 

The V 2 results are similar between AMPT with Debye 
cr = 3 mb and MPC with Debye a = 5.5 mb, both cor¬ 
responding to a similar opacity (A^coii) (see Table |T]). In 
Fig. I^a) we show instead the MPC results with isotropic 
tj = 5.5 mb to compare to the AMPT results. The results 
are qualitatively similar while the freezeout (^ 2 ) is larger 
in MPC. MPC results with a larger opacity are shown in 
Fig. I^b). The escape picture still holds at large opaci¬ 
ties. However, the probabilities to escape at small A^coii 
are now smaller, so the escape contribution to the final 
overall V 2 is smaller, and the remaining active partons 
have mostly positive V 2 - 

Figure]^ shows the approximate average transverse ra¬ 
dial velocity of partons in AMPT, /3j_ = (fj_ • p±) where 
f± and p± are the transverse radial position and momen¬ 
tum unit vectors, as a function of Ncoii- The /3± of all 
partons (thick solid curve) at A^coii = 0 is not exactly 
zero because partons can form only after a finite forma¬ 
tion time over which a parton’s displacement depends 
on its momentum. The freezeout partons (thick dashed 
curve) at A^coii = 0 have a large /3j_. This strong space- 
momentum correlation is due to the anisotropic escape 
mechanism, but different from a collectivity that repre¬ 
sents a common collective flow velocity achieved only 
via interactions. On the other hand, there are space- 
momentum correlations for all partons (thick solid curve) 
at any given A^coii > 0 and the correlation increases with 
Afcoii; these correlations are good indicators of collective 
flow. Some of these partons freeze out at a given A^coii; 
the additional (3± for these freezeout partons is the effect 
of the anisotropic escape mechanism. 

One question is whether the hydrodynamic-type col¬ 
lective flow of the partons is important for the final V 2 - 



FIG. 4: Parton /3u = ' P±) as a function of Aeon in 

Au-|-Au collisions. Both normal (thick curves) and azimuth- 
randomized (thin curves) AMPT results are shown. The solid 
curves are those for all partons after suffering Aooii collisions, 
and the dashed curves are for freezeout partons (partons that 
freeze out after suffering exactly Acoii collisions). 


Thus we did test calculations with no collective flow by 
randomizing the outgoing parton azimuthal directions af¬ 
ter each parton-parton scattering. The system continues 
to evolve in AMPT, but the evolution is different from 
the original one. The /3 _l from this modified evolution is 
shown in Fig.[^ The all-parton fl± is now zero because of 
the randomization, and the freezeout parton fl± is non¬ 
zero purely due to the anisotropic escape mechanism. 

We show in Fig. the V 2 of all partons and freeze¬ 
out partons from this azimuth-randomized AMPT by the 
thin solid and dashed curves for Au-|-Au and d-|-Au col¬ 
lisions. In the randomized case, the parton azimuthal 
angles are randomized and hence their V 2 is zero; thus 
the final-state freezeout anisotropy is entirely due to the 
anisotropic escape mechanism. For comparison, the V 2 
results from the normal AMPT (already shown in Fig. 
for Au-|-Au) are superimposed in Fig. as the thick solid 
and dashed curves, where the all-parton V 2 is slightly 
positive and the freezeout parton V 2 is much higher. The 
gain in V 2 by the freezeout partons is due to the escape 
mechanism. The gain in the normal AMPT results is 
slightly different from that in the azimuth-randomized 
results. This is not surprising because the anisotropies 
in the escape probability differ in these two cases: in the 
former case the parton p ±’s are correlated with their f ±’s 
while in the latter case the parton p ±’s are random. 

We have shown mostly results for 200 GeV Au-|-Au 
collisions at medium impact parameters. There seems 
to be no qualitative difference between the behaviors 
in Au-I-Au and d-|-Au collisions (c.f. Fig. [^. Although 
we focused on V 2 ^ the same qualitative conclusions hold 
for V 3 as well-see Fig. where V 3 is shown similar to 
Fig.[5];a). This suggests that the development mechanism 
of anisotropies in transport models is universal. We note 
that a fixed Acoii value does not correspond to partons 
at identical time but rather a convolution over time. We 
have also studied results as a function of time instead of 
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FIG. 5: Parton V 2 as a function of A^coii in (a) Au+Au and 
(b) d+Au collisions. Both normal (thick curves) and azimuth- 
randomized (thin curves) AMPT results are shown. The solid 
curves are for all partons after suffering Acoii collisions, and 
the dashed curves are for freezeout partons. 


FIG. 7: Gumulative V 2 of all partons (see text) as a function 
of A^coii in (a) Au-|-Au and (b) d-|-Au collisions. Both normal 
(thick curves) and azimuth-randomized (thin curves) AMPT 
results are shown. 



FIG. 6: Parton vs as a function of Acoii in Au-|-Au colli¬ 
sions. Both normal (thick curves) and azimuth-randomized 
(thin curves) AMPT results are shown. The solid curves are 
for all partons after suffering Aeon collisions, and the dashed 
curves are for freezeout partons. 


Acoii, and our qualitative conclusions remain unchanged. 

Experimentally, only the final-state anisotropy is mea¬ 
sured, integrated over all evolution time. It is thus in¬ 
teresting to examine the cumulative V 2 of all partons up 
to Acoii collisions, including those that have frozen out 
with iVcoii or fewer collisions and those that will suffer 
further collisions, as a function of iVcoii. This is shown 
in Fig. The cumulative V 2 increases with JV^oii and 
starts to saturate after approximately 10 and 3 collisions 
in mid-central Au-|-Au and central d-|-Au collisions, re¬ 
spectively. The asymptotic V 2 value at Acoii —t oo would 
be the final-state parton average (^ 2 ) in the events. 

Table|T]lists the (Acoii) and (^ 2 ) of all final partons from 
normal and azimuth-randomized results by both AMPT 
(varying with Ug kept fixed) and MPC (isotropic 
a). The {V 2 ) with Debye u of 3 mb and 5.5 mb from 
AMPT and MPC, respectively-so that they have simi¬ 
lar opacity-are consistent with each other as mentioned 


previously. It is known that V 2 is larger for an isotropic 
da/dt than for a Debye screened one (at the same a) be¬ 
cause the latter is more forward-backward peaked, but 
the (r' 2 )Rndm/('i' 2 ) ratios as well as the opacities are al¬ 
most the same. It should be noted that there may be 
issues with causality in AMPT for large a. Also note 
that (Ncoii) is larger in the randomized case because the 
randomization tends to destroy the preferred outward di¬ 
rection of partons. 

Discussions. The unique finding of our study is that 
partons in transport models escape (freeze out) from the 
collision zone with positive U 2 , even those partons that do 
not interact at all, mainly due to the anisotropic escape 
probability. This escape mechanism contributes to the 
majority of the final V 2 at small to modest opacity. The 
traditional picture of low p_L particles accumulating V 2 
after multiple collisions seems to play only a minor role. 
High-p_L anisotropy is generally believed to result mostly 
from the escape mechanism [531 HI] > which we have also 
verified within our work. Our study indicates that the 
escape mechanism is at work at both high and low p_l; 
there appears to be no fundamental difference in the V 2 
development of high- and low-pj^ partons. Thepj^ depen¬ 
dence of r! 2 -also captured in AMPT-is less obvious from 
the escape mechanism. It is partly borne out of the fact 
that higher-p_L partons freeze out after fewer collisions 
(c.f. Fig. [^a)), and hence possess larger V 2 (c.f. Fig. |^. 
The LDL calculation of Ref. |55| found that the central¬ 
ity dependence and the magnitude ofp_L-averaged elliptic 
flow were well described, but the shape of the p± depen¬ 
dence was siginificantly off at RHIC. 

It has generally been thought that the reason why 
AMPT describes the bulk experimental data well m 
is because AMPT has large opacity and thus has ap¬ 
proached hydrodynamics. Our study indicates that this 
interpretation is incorrect and in fact the opacity in 
AMPT is only small to modest. If one takes the (u 2 )Rndm 
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TABLE I: (A^coii) and {V 2 } of all final partons from normal (first number in each column) and azimuth-randomized AMPT 
(Debye a) and MPC (isotropic a) results (second number). The d-l-Au impact parameter is 6 = 0 fm. 



AMPT d-l-Au 

AMPT Au-|-Au {b = 7.3 fm) 

MPC Au-tAu {b = 8 fm) 

(J 

3 mb 

3 mb 

20 mb 

40 mb 

60 mb 

5.5 mb 

20 mb 

40 mb 

60 mb 

(Neon) 

1.2 1.4 

4.6 5.8 

13 22 

17 32 

20 39 

4.7 5.4 

17 23 

35 52 

53 83 

{V2) 

2.7% 2.5% 

3.9% 2.7% 

5.9% 2.7% 

6.0% 2.3% 

5.7% 2.0% 

5.5% 3.5% 

8.6% 3.5% 

9.8% 3.0% 

10% 2.6% 

(r'2)Rndm/ {V 2 ) 

93% 

69% 

46% 

38% 

34% 

64% 

41% 

31% 

26% 


from the azimuth-randomized results as estimate of the 
escape contribution to the {V 2 ) from the normal results, 
then its contribution in semi-central Au-|-Au collisions 
is ^70% with modest opacity (Acon)=4-5. As opacity 
increases, the hydrodynamic-type collective flow contri¬ 
bution increases, but rather slowly. The system is still 
far away from asymptotic hydrodynamic behavior even 
with the unrealistically large opacities we have studied. 
It is found that this picture is qualitatively the same for 
V3- 

The space-momentum correlations of freezeout partons 
in the transport models are largely due to the escape 
mechanism. This results in “surface emission,” [24] where 
a parton freezes out depending on its momentum and 
position, which determine its escape probability at that 
point of evolution. It thus creates a space-momentum 
correlation even in the absence of hydrodynamic collec¬ 
tive flow, as demonstrated by the azimuth-randomized 
AMPT results. It is important here to distinguish be¬ 
tween space-momentum correlation and hydrodynamic 
collective flow, where the latter means a collective motion 
that is generated by interactions so that particles convert 
part of their energy into their common motion (e.g. parti¬ 
cles in nearly local thermal equilibrium moving on top of 
a common velocity field). There is indeed a finite collec¬ 
tive flow in AMPT (the thick solid curve in Fig.|^. This 
radial flow is presumably generated by hydrodynamic- 
type interactions and pressure gradient. The azimuthal 
modulation of this radial flow is the anisotropic flow Vn of 
partons (the solid thick curves in Figs. mm and . It is 
the Vn of these active partons that is the most relevant for 
the standard hydrodynamic flow description, or the col¬ 
lective properties of the sQGP. The radial flow may be a 
viable discriminator: collectivity generates extra p± but 
the escape mechanism does not. Experimentally, how¬ 
ever, radial flow may be difficult to measure as it requires 
precise knowledge of initial-state transverse energy pro¬ 
duction. The centrality dependence of Vn may serve as 
another possible discriminator, since hydrodynamics and 
LDL give different predictions |2T1 [22| . 

Hydrodynamics has been successfully applied to heavy 
ion collisions |2], and a small viscosity to entropy den¬ 
sity ratio (rj/s), close to the conjectured quantum limit 
of l/47r, has been extracted [25]. Hydrodynamic evolu¬ 
tion is typically stopped in a calculation when the lo¬ 
cal energy density or temperature reaches a given value. 


Particle production is then modeled by the Cooper-Frye 
formalism |26j . The escape mechanism, on the other 
hand, is driven by the chance of no further interaction 
^escape = exp (— J padtj, i.c., a non-local quantity that 
involves the entire future density (p) evolution of the sys¬ 
tem, and is not obviously captured by the Cooper-Frye 
prescription. If hydrodynamically driven collective flow 
is indeed a small contribution to the experimentally mea¬ 
sured anisotropy, then the extracted rj/s is severely un¬ 
derestimated. Since escape is inevitable for our transient 
collision systems, it is imperative to examine the possi¬ 
ble role of the escape mechanism in the hydrodynamics 
framework. Previous studies have shown that continu¬ 
ous particle emission instead of sudden freezeout in hy¬ 
drodynamics can have important implications for pion 
interferometry |27j . 

In summary, we have studied the development of az¬ 
imuthal anisotropy 'u„(n = 2,3) in AMPT and MPC as 
a function of the number of collisions A^coii that a par- 
ton suffers in Au-I-Au and d-|-Au collisions at = 

200 CeV. It is found that the majority of comes from 
the anisotropic escape probability of partons, and this 
picture applies similarly to partons at both high p± and 
low p±. The anisotropic flow of partons as a result of 
parton-parton interactions or hydrodynamic-type pres¬ 
sure gradient is found to be small in transport models. 
This part of anisotropy becomes more important with 
increasing cross-section a, but the system in transport 
models is still far from the asymptotic hydrodynamic be¬ 
havior even with the unrealistically large cross-sections 
that we studied. The escape mechanism is dominant 
because the evolution in transport models is relatively 
dilute ((A^coii)=4-5 at ct ^ 3 mb). This is distinctly dif¬ 
ferent from evolution near the hydrodynamic limit, where 
collectivity is generated by a large number of collisions. 
Such hydrodynamic-type collectivity and anisotropic flow 
is what one would regard as the cornerstone of the sQCP 
paradigm. Our results, however, suggest that the as¬ 
sumption of hydrodynamics is not imperative. 
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